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ABSTRACT 

Aims. We investigate the effect of ram pressure stripping (RPS) on several simulations of merging pairs of gas-rich spiral galaxies. 
We are concerned with the changes in stripping efficiency and the time evolution of the star formation rate. Our goal is to provide an 
estimate of the combined effect of merging and RPS compared to the influence of the individual processes. 

Methods. We make use of the combined N-body/hydrodynamic code GADGET-2. The code features a threshold-based statistical 
recipe for star formation, as well as radiative cooling and modeling of galactic winds. In our simulations, we vary mass ratios between 
1:4 and 1:8 in a binary merger. We sample different geometric configurations of the merging systems (edge-on and face-on mergers, 
different impact parameters). Furthermore, we vary the properties of the intracluster medium (ICM) in rough steps: The speed of 
the merging system relative to the ICM between 500 and 1 000 km s -1 , the ICM density between 10~ 29 and 10 -27 g cm -3 , and the 
ICM direction relative to the mergers’ orbital plane. Ram pressure is kept constant within a simulation time period, as is the ICM 
temperature of 10 7 K. Each simulation in the ICM is compared to simulations of the merger in vacuum and the non-merging galaxies 
with acting ram pressure. 

Results. Averaged over the simulation time (1 Gyr) the merging pairs show a negligible 5% enhancement in SFR, when compared 
to single galaxies under the same environmental conditions. The SFRs peak at the time of the galaxies first fly-through. There, our 
simulations show SFRs of up to 20 M 0 yr -1 (compared to 3 M 0 yr~' of the non-merging galaxies in vacuum). In the most extreme 
case, this constitutes a short-term (< 50 Myr) SFR increase of 50 % over the non-merging galaxies experiencing ram pressure. The 
wake of merging galaxies in the ICM typically has a third to half the star mass seen in the non-merging galaxies and 5% to 10 % less 
gas mass. The joint effect of RPS and merging, according to our simulations, is not significantly different from pure ram pressure 
effects. 

Key words. Galaxies:clusters: general - Galaxies:ISM - Galaxies:interaction - Galaxies:spiral - Galaxies: clusters:intracluster medium 
- Methods: numerical 


1. Introduction 


The changes galaxies undergo during their lifetimes are caused 
by complicated interplays between many physical processes 
(Lackner & Gunn 2013). Isolating just a single contributor to 
galaxy evolution as the dominant one might not be possible (Mi- 
hos 2003), and so the focus in theoretical work is generally shift¬ 
ing toward a more wholesome description of galaxy evolution. 
However, two physical processes, ram pressure stripping (RPS) 
and galaxy interaction, are thought to be the potent driving forces 
behind galaxy evolution. 

Ever since Gunn & Gott (1972) introduced the idea of RPS 
of studying galaxy evolution, this effect has received more and 
more attention: from early observations (e.g., Irwin et al. 1987; 
Cayatte et al. 1990) to distant clusters and the existence of star 
formation in the stripped wake (e.g.. Sun et al. 2007; Hester et al. 
2010) that Abadi et al. (1999) and Farouki & Shapiro (1980) 
pioneered in simulating RPS numerically, essentially confirm¬ 
ing and refining the theoretical estimates of Gunn & Gott. In 
their simulations, Vollmer et al. (2001) and Schulz & Struck 
(2001) found a mechanism for enhancing star formation in the 
inner region of galaxies not affected by RPS and discussed the 
role of RPS in the Butcher-Oemler effect (e.g.. Butcher & Oem- 
ler 1978, 1984). Simulations have been carried out using La- 
grangian methods, for example, in order to quantify the RPS 
effects on internal kinematics (Rronberger et al. 2006), star for¬ 


mation (Kronberger et al. 2008, Kapferer et al. 2009 hereafter 
K9), and stellar bulge size (Steinhauser et al. 2012 hereafter SI). 
Among others, Roediger & Briiggen (2006, 2007) used Eulerian 
grid codes, which complement the Lagrangian schemes. Ton- 
nesen & Bryan (2012) and later Roediger et al. (2014) explicitly 
target star formation in RPS tails. 

Similarly, galaxy mergers are known to incite star formation 
activity (e.g., Elbaz & Cesarsky 2003; Joseph et al. 1984). The 
remarkable early simulations of Toomre & Toomre (1972) con¬ 
nected merger activity to the formation of tidal bridges and other 
features, leading to redistribution of gas to the core and thus star 
formation. The simulations of Noguchi (1988) show that merg¬ 
ers could indeed trigger such an inflow. This general theory of 
merger-induced star formation has since been supported by a 
myriad of publications. More recently, Hopkins et al. (2011); 
Bournaud et al. (2010) have focused on refining these models in 
an attempt to improve the resolution of individual giant molecu¬ 
lar clouds as clumps of star formation. Additional effort has been 
put into modeling star formation processes correctly, starting 
with the Kennicutt-Schmidt law (see Schmidt 1959; Kennicutt 
1998). Various improvements on these semi-empirical laws have 
been made since Schmidt’s publication (see, e.g., Elmegreen 
2011 for a general review). 

Owing to the high velocity dispersion in massive clusters, 
it has been long thought that only high-speed fly throughs and 


Article number, page 1 of 12 



Table 1: Important parameters used in all simulations (see 
Springel & Hernquist 2003; Springel 2005). 


Parameter 

Symbol 

Value 

Supernova fraction 

P 

0.1 

Evaporation factor 

A 

1000 

Supernova temperature 

7sn 

10 8 K 

Cold cloud temperature 

T c 

1 000 K 

Threshold density 

Pth 

2.89 x 10~ 25 g cm -3 

Star formation timescale 

t* 

*0 

2.1 Gyr 

Wind efficiency 

n 

2 

Wind energy fraction 

X 

0.5 


fly bys would occur. Thus mergers, requiring low relative ve¬ 
locities between the merging partners, seemed very rare (e.g., 
Ostriker 1980). However, mergers do occur, because clusters 
form hierarchically and merge with infalling groups of lower 
dispersion (Mihos 2003). Gnedin (2003) simulated that on av¬ 
erage each galaxy has about one encounter with relative veloc¬ 
ity < 500 km s -1 . Depending on the criterion for cluster mem¬ 
bership, the fraction of cluster galaxies that originate in groups 
varies (De Lucia et al. 2012). McGee et al. (2009) estimate that, 
depending on cluster mass, between 30% and 45% of galaxies 
(M* > 10 9 M 0 h” 1 ) are accreted from groups that are more mas¬ 
sive than 10 13 M 0 h _1 . Berrier et al. (2009) give a fraction of 
10% and find that ~30% of cluster galaxies are accreted with at 
least one luminous companion. 

All this motivates a closer look at mergers in a cluster en¬ 
vironment. To date, most theoretical work has regarded galaxy 
interactions and RPS independently. This approach is viable for 
disentangling the different phenomena, and it leads to an under¬ 
standing of each phenomenon individually. In real galaxy clus¬ 
ters both processes occur at the same time. It cannot be said in 
advance whether both processes working together have syner¬ 
getic effects, notably on star formation and stripping fraction. In 
this paper we expand the work of Kapferer et al. (2008), here¬ 
after K8, to a larger set of simulations and refine the techniques 
used there. That means that we explore a grid of merger and RPS 
parameters and make cross-comparisons to scenarios where one 
of the phenomena (merger or ram pressure) is missing. 


2. Simulation 

We use the smoothed-particle hydrodynamics (SPH) and N-body 
code GADGET-2, developed by Volker Springel, for all simu¬ 
lations in this work. Detailed information on the code can be 
found in Springel (2005). The N-body part responsible for the 
gravity calculation uses a geometric Barnes & Hut (1986) oct- 
tree grouping with monopole moments in the force expansion 
and an error-limiting cell-opening criterion. Long-range forces 
are calculated according to the TreePM method (e.g.. Xu (1995)) 
using a Fourier-based approach. This is especially efficient for 
the periodic boundaries of our ICM setup (described below). 
GADGET-2 uses an energy- and entropy-conserving formula¬ 
tion of SPH (Gingold & Monaghan 1977; Lucy 1977) as pre¬ 
sented in Springel & Hernquist (2002). 

Implemented physics include radiative cooling (Katz et al. 
1996) and, as described in Springel & Hernquist (2003), star 
formation, stellar feedback, and galactic winds. Relevant pa¬ 
rameters of our simulations are summarized in Table 1. Most 
notably Table 1 includes parameters for the statistical star forma¬ 
tion recipe. Star formation is implemented as a subgrid model, 
since single stars cannot be resolved individually. Star formation 


sets in when the density of a gas particle exceeds a threshold p t h. 
The density of gas particles in star-forming regime is partitioned 
as a cold and a hot part, p gas = p c + Ph- Also considering the 
newly formed stars, the masses of these two phases evolve ac¬ 
cording to 


d Pc 



SF 



+ 


Anet(p/i> U h ) 

Uh ll c 


d __ _ n Pc , \ 0 P C 1 f k t „ \ 

Al _Ph — ‘ _ 'XietCP/;, Uh) 


dt 


r//j n c 


SN heating 


rad. cooling 


New stars are formed out of cold gas on a characteristic 
timescale t*. It is assumed that massive stars (M > M 0 , con¬ 
stituting a fraction /3 = 0.1 in a Salpeter IMF) immediately blow 
up as supernovae with their gas added to the hot phase. Fur¬ 
thermore, cold gas is evaporated thanks to star formation with 
efficiency A and radiative cooling transfers hot gas to the cold 
phase. In a similar way, energy exchange between the gas phases 
is modeled. The implementation of this star formation recipe is, 
however, done in a statistical fashion to avoid the calculation of 
the mass exchange processes. In the case of self-regulated star 
formation, the hot gas phase evolves toward an equilibrium tem¬ 
perature. This also determines the fraction in cold clouds and 
the star formation rate can be calculated as M* = (1 — j8)^“- 
Consequently, new star particles are spawned from a given gas 
particle with a probability proportional to 1 - exp in a 

time step of length At. These particles conserve the phase space 
properties of the progenitor gas particle. For more details, see 
Springel & Hernquist (2003). 

In our simulations, we vary geometric properties and the 
mass ratio of the galaxy mergers, ICM densities and relative ve¬ 
locities to investigate the combined effects of RPS and minor 
mergers on SFR and stripping. Overall, 20 core simulations with 
different combinations of parameters were conducted, capturing 
the first two encounters of the constituent galaxies and spanning 
a simulation time of 1 Gyr each. Including simulations done 
for comparison purposes, the total number of different individual 
runs is 58. Table 3 summarizes our initial conditions. To have an 
acceptable number of simulations, we do not sample the whole 
grid of parameters. Rather we vary a standard case referred to 
as wstd in only one parameter at a time. The parameters used 
in Table 3 are described in the following subchapters along with 
the galaxy and ICM model. 

For the naming convention in Table 3, the different simula¬ 
tions are named after their type and the main parameter(s) that 
deviate from the standard case labeled wstd. For instance, the 
name “wp27” is used for a simulation of an ICM density of 
Picm = 10 27 g cm -3 . Otherwise, wp27 is identical to the stan¬ 
dard case in its parameters. Simulations starting with the letter 
“w” are conducted in our “windtunnel” setup described in Sect. 
2.2. These simulations are meant to model merging galaxies in 
ICM wind. The starting letter “i” indicates single (non-merging) 
galaxies in the ICM, which are necessary for later comparison. 
Mergers in vacuum have an “m” in the beginning of their desig¬ 
nations. A simulation starting with none of these letters refers to 
a single galaxy in vacuum. 


2.1. Galaxy model 

The galaxies are created by the initial conditions generator used 
in Springel et al. (2005), kindly provided by Volker Springel. 
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Table 2: Initial galaxy properties for the largest galaxy in our 
simulations. Those of the smaller galaxies are printed in brack¬ 
ets (where different). Values are before evolution of 1 Gyr, i.e., 
before self regulation. The first seven are free parameters. 


Galaxy property symbol value (smaller galaxies) 


Halo concentr. 

c 

7 


Circular velocity 

V 200 

180 km s'" 1 

(113, 99, 90) 

Spin parameter 

A 

0.06 


Spin fraction 

jd 

0.02 


Disk scale height 

hd 

0.1 



# gas particles 

N g 

300000 

(75 000, 50000, 37 500) 

# disk particles 

N d 

-N 

3 iy 8 





# halo particles 

N h 

-N 

3 iy 8 





Disk scale length 

l 

6.52 

kpc 


(4.11,3.59, 2.90) 

Gas mass 

Mg 

4.06 

X 

10 9 

M© 

(1.01,0.68,0.51) 

Stellar mass 

M d 

2.71 

X 

10 10 

M© 

(0.68, 0.45, 0.34) 

Total mass 

Tftot 

1.75 

X 

10 12 

M© 

(0.44, 0.29, 0.21) 

Gas resolution 


1.36 

X 

10 4 

1^0 

particle 


Stellar resolution 


1.36 

X 

10 5 

Mq 

particle 


Halo resolution 


4.32 

X 

10 6 

M© 

particle 



The theoretical groundwork for the generator was described in 
Mo et al. (1998). The generator uses a radial Hernquist (1990) 
profile of scale length a for the dark matter spatial density: 

Pdm(r) — — - — -3 

In r(r + ay 


It relates to the more common NFW (Navarro et al. 1996) profile 
via a - r s V2[ln(l + c) - c/(l + c), where c = r 2 oo/r s is the con¬ 
centration parameter with r s the scale length of the NFW halo, 
and r 2 oo the radius enclosing a mean dark matter density of 200 
times the critical density. Gas and stars in the disk follow an 
exponential surface density profile: 


Z*, g (R) 


M*lg „-R/l 
2nP 


The profile falls with disk scale length I. This in turn is computed 
from the angular momentum as a function of the spin parameter 
A, circular velocity V 200 , and halo concentration c. In the case 
of gas particles, the vertical mass distribution is sampled from 
the overall potential using the equation of state. In the case of 
star particles, it is modeled as an isothermal disk of relative disk 
height hj. Additional parameters are gas and stellar disk mass 
fractions, nij and m g , relative to total mass and disk mass, re¬ 
spectively. The parameters of the galaxies used here are summa¬ 
rized in Table 2. The parameters in Table 2 were mainly selected 
for their stability. In 4 Gyr of evolution, they showed no morpho¬ 
logical peculiarities and yielded no unexpected results in terms 
of SFR. We used pure disk galaxies, i.e., without a bulge com¬ 
ponent or Hl-disk. The galaxies in our simulations differ only 
in their V 200 parameter, representing the circular velocity of the 
halo at r - t^oo- We varied only the mass of the smaller merging 
partner as seen in Table 2. This gives mass ratios of p = i, p — ~ 
and p = | compared to the largest galaxy in use. In all, the cho¬ 
sen parameters are fit to represent gas-rich spiral galaxies. The 
baryonic mass of our largest galaxy constitutes 2% of its total 
mass, where 15% is gas. The particle numbers of the different 
components are chosen empirically by the same method as used 
by Kapferer et al. (2005). There it was found that star formation 
and gas distribution are essentially independent of the gas mass 



(a) Initial gaseous disk (b) Gaseous disk after 1 Gyr 

Fig. 1: Particle distribution in disk galaxies created with the gen¬ 
erator of Springel et al. (2005). After 1 Gyr of isolated evolution 
spiral structure has developed and star formation self regulated 
(not shown). 


resolution, as long as it is below 10 6 pa W. le . The collisionless disk 
and halo particles can be sampled at a higher mass per particle 
because they are numerically less problematic. 

The initial-condition generator provides a density and veloc¬ 
ity distribution of particles. We evolve the created galaxies for 
1 Gyr in isolation. In this time, a spiral arm structure forms (see 
Fig. 1) and 4.8 x 10 8 M© of stars are produced in the largest 
galaxy. The radial density structure of our highest mass galaxy 
is very similar to the one described in SI (their Fig.l), but fur¬ 
ther extended owing to the higher circular velocity in our simu¬ 
lations. After isolated evolution, the galaxies are paired up into 
binary mergers and, in case of the main runs, inserted into the 
ICM environment. 


2.2. ICM setup 

GADGET-2 features only cubical simulation domains with pe¬ 
riodic boundaries. For typical galaxy velocities within a cluster 
(1 000 km s 1 ), this poses a computational problem: To capture 
at least the first two encounters of a merging pair of galaxies, 
the domain needs a size of about L — 1 Mpc cubed. To simu¬ 
late the ICM with GADGET-2, one could fill this domain with 
gas particles of sufficient mass resolution. To avoid numerical 
artifacts caused by the SPH scheme, the mass resolution of ICM 
particles has to match that of the galaxie’s gas (Ott & Schnetter 
2003). Accurate simulations of both collisionless dynamics and 
smoothed particle hydrodynamics require a certain mass resolu¬ 
tion. as does the star formation recipe. To fill the whole domain 
with particles of this mass resolution is computationally too ex¬ 
pensive. We adapted the approach in K9 with a layered setup 
with resolution decreasing toward the outer layers and periodic 
boundaries. 

Our setup consists of an inner elongated “windtunnel” of 
high-resolution particles encompassed by particles with lower 
resolved masses. These coarser resolved particles keep the in¬ 
ner region stable. To prevent the setup from dissipating, we use 
periodic boundary conditions. In contrast to K9, our setup intro¬ 
duces an additional layer, a sleeve region around the innermost 
part, which drops slightly in resolution. This improves stability 
of the inner region and allows for a steeper resolution drop in 
the large outer region, thus saving computation time. The setup 
is shown in Fig. 2. The overall ICM density does not change 
by more than 2 % in the inner region of our simulations. Dissi¬ 
pation and density fluctuations are suppressed better than in K9. 
The mean distance of particles that have escaped the inner re¬ 
gion from the region’s edges is up to 20% less in our setup than 
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Fig. 2: Three-layered ICM setup. Shown is the region of high 
resolution (inner border, hxhxL) in which the simulated merger 
takes place, the encompassing sleeve region of lower resolution 
(border S x S x L) and the rest of the simulation domain (out¬ 
ermost cube of edge length L). The particle’s number density 
is indicated in shades of green. Galaxies move along the arrow 
inside the region of high-mass resolution. 


with K9, while their number decreases by 30%. Depending on 
the merger geometry, the “windtunnel” typically has a size of 
120 x 120 x 1 000 kpc with a sleeve width of 60 kpc. 1 This 
gives a time of 1 Gyr for the galaxies to pass through. 

The high-resolution region is large enough to hold all of the 
galaxies’ gas and most of the wake at all times during the merger. 
The cube size is sufficient for the whole extent of the galaxies’ 
dark matter halo. It must be noted, however, that the periodic 
images of the galaxies still contribute to the overall potential. 
These are in principle unphysical and introduce small force er¬ 
rors. The ICM throughout this work is modeled with a tem¬ 
perature of 10 7 K, which is consistent with typical ambient gas 
temperatures in a cluster. 

In other ram pressure simulations, SI and Jachym et al. 
(2007, 2009) find little dependence of SFR and stripping frac¬ 
tion on the direction of ICM flow with respect to disk orienta¬ 
tion. Nearly edge-on cases (within « 20°) only had significant 
deviations in stripping fraction from the cases where galaxies en¬ 
countered the ICM with their disks oriented face-on to the ICM. 
However, the differences in SFR were negligible in all cases. 
Based on these results, and owing to other restrictions discussed 
in Sec. 2.3, we chose to sample mostly face-on cases, with two 
edge-on ICM scenarios for completeness. The ICM densities in 
our simulations range from 10 27 to 10 -29 g cm -3 . Adopting a 
standard /j model 2 * * for the ICM density, these values correspond 
to cluster-centric distances of 0.5, 0.8, 1.8, 2.6, and 5.6 Mpc 
for Coma, putting our standard case at 0.7 r^° ma . As such they 
cover typical densities from the outskirts to the regions where 
ram pressure dominates our simulations. For Virgo our standard 
case would correspond to a distance of about 1 Mpc from the 


1 Corresponds to h x h x L in Fig. 2. Sleeve width: = 60 kpc. 

2 /3-Model of Cavaliere & Fusco-Femiano (1978): PicmM = Poll + 

(r/r c ) 2 ]' 3/20 . Parameters for the Coma cluster (Mohr et al. 1999) are: 

/3 = 0.7, r c = 0.26 Mpc andpo = 6.3 x 1CT 27 g cm -3 . 



Fig. 3: Geometric merger parameters as used in our simula¬ 
tions.With d, the inital separation, b the inital displacement, as 
well as p and y the inclination of the larger and smaller galaxies 
with respect to the orbital plane. 


center and 0.3 Mpc for the highest density in our simulations 
(Schindler et al. 1999, Fig. 1 lb). 

The strength of ram pressure, as discussed by Gunn & Gott 
(1972), is given by 

Fram = PlCM 

where P lam is the pressure due to the relative movement of galax¬ 
ies with respect to the ICM at speed v re i- A gas density of the 
ICM, picm = const., is assumed in our simulations. The amount 
of stripped gas in this approximation is then determined by the 
radius at equilibrium between ram pressure and restoring gravi¬ 
tational force. A similar estimate for the stripping timescale can 
be made by comparing acceleration due to ram pressure with 
restoring acceleration due to gravity. 

It was found, e.g., by K9 that SFR has a stronger dependence 
on picm than on v re i. Consequently, we sampled five different 
ICM densities in contrast to only two speeds within the ICM 
(see Table 3). 


2.3. Merger setup 

For the most general case of a binary merger, the centers of mass 
of the two galaxies approach each other at approximately Kep- 
lerian orbits. The angular momenta of their disks with respect to 
the orbital plane can have arbitrary orientations. The ICM flow 
may come from an arbitrary direction as well. 

However, the treatment of ICM in our simulations limits the 
range of some of the individual geometric parameters of the 
merger; that is, we can only simulate those geometries that fit 
within the high-resolution region of the ICM shown in Fig. 2. 
For convenience, we align the orbital plane with our “windtun¬ 
nel”. To capture the whole merging process from the start, the 
initial separation of both disks, d, must be at least 100 kpc. The 
number is supported by both our own simulations and observa¬ 
tions (Ellison et al. 2010; Li et al. 2008; Lambas et al. 2003). 
At this distance, gravitational interaction can be neglected for 
star formation and stripping, and the whole merging process is 
captured. 

Lambas et al. (2003) found significant enhancement in star 
formation only in galaxies at relative velocities v < 100 km s _1 , 
which is the highest relative initial velocity between two galax¬ 
ies, which we chose to adapt in our simulations. The effective 
set of geometrical parameters are illustrated in Fig. 3. 
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Table 3: Designation and properties of the different simulations conducted in this work. 


Parameter 

(where different) 

wstdirf 

wp4(r) a 

wp8(r) a 

wp21 

wp528 

wp529 

Simulation designation 

wp29 wvr500 wvm50 

wvm!5 

wblti 

wb 50 

wy45 

wy 90 

wee 

w/3y° 

Unit 


1:6 

1:4 

1:8 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

1 

n d 

PlCM 

1(T 28 

- 

- 

10~ 27 

5 x 10~ 28 

5 x 10- 29 

10“ 29 

- 

- 

- 

- 

- 

- 

- 

- 

- 

g ctrT 3 

Vrel 6 

1000 

- 

- 

- 

- 

- 


500 

- 

- 

- 

- 

- 

- 

- 

- 

km s’ 1 

Vmer f 

100 

- 

- 

- 

- 

- 

- 

- 

50 

75 

- 

- 

- 

- 

- 

- 

km s” 1 

b% 

30 

- 

- 

- 

- 

- 

- 

- 

- 


10 

50 

- 

- 

- 

- 

kpc 

7 h 

0 (180) 

0 (180) 

0 (180) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

45 

90 

90 

90 b 

° (deg.) 

/? h 

0 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

90 

90 b 

° (deg.) 

d' 

100 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

kpc 

merger 
large galaxy 
small galaxy 1 

mp6 

ipl 

i/j6 

m/.i4 

ip4 

mp8 

i/x 8 

ip\p21 

ip6p21 

;plp528 

;p6p528 

;plp529 

;p6p529 

ip 1 p29 
;p6p29 

ip 1 vt'500 
/p6v;500 

mvm50 

mvml5 

mblO 

mb50 

my45 

ip6y45 

my90 

ip6y90 

mee 

in 6y90 

on ►r 3 

^ 7=2 2 

SO & 

O O 

n.A. 

n.A. 

n.A. 


Notes. (a) Retrograde partner simulation, i.e., with y = 180° end with “r”, e.g. in wp4r the two disk are rotating in opposite directions (b) The galaxies in wfiy merge edge-on, both with ICM face-on 
(c) p, mass ratio between the two galaxies (d ' p K M, density of the ICM in the simulation domain (e) v re i, center of mass velocity of the merging galaxies relative to the ICM ® v mer , initial velocity 
of the galaxies relative to each other (8) b, impact parameter of the merging galaxies <hl /i and y, inclination angles to the orbital plane of the merger’s more resp. less massive galaxy (l) d, initial 
separation of the galaxies (1) In the last three rows, the merger simulation in vacuum (starting letter “m”) and the two single-galaxy simulations in ICM (“i”) are listed. 
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Fig. 4: Overall SFR of the reference simulation wstd and its 
comparison simulations, ip 6 + i/j.1 and mp 6. Shown is the added 
SFR of the single galaxies in ICM and in vacuum (upper and 
lower dashed lines, respectively), as well as the SFR for the same 
galaxies merging in vacuum and in an ICM environment (solid 
lines). The strong peak in vacuum and ICM merger SFRs coin¬ 
cides with the galaxies first encounter at t — 0.4 Gyr. 

3. Results 

3.1. Evolution of star formation rate 

In this section we investigate the impact of the simulated ram 
pressure wind on our minor merger’s star formation rate (SFR). 
Figure 4 shows the evolution of SFR over 1 Gyr of simulation 
time for our reference case, wstd, and its three comparison runs 
(non-merging galaxies in vacuum and ICM, merger in vacuum). 
Similar plots are provided for the most interesting cases simu¬ 
lated. In the plots, we show the total SFR, including stars form¬ 
ing in the ram-pressure-induced wake behind the galaxies, as 
well as those in the galactic disks. In the next section (3.2), we 
discuss the distribution of newly formed stars between the two. 

In most plots, the SFR shows two distinct peaks at about 25 
and 400 Myr into the merger’s evolution. The initial peak in the 
SFR plots (at ~ 25 Myr) should be disregarded. It occurs when 
the galaxies first encounter the ICM of our setup. In reality they 
would gradually fall into a cluster’s ICM as opposed to a sud¬ 
den encounter, when inserted into the “windtunnel” (Sec. 2.2). 
The disturbance caused by that, however, is overcome after the 
first 80 Myr of evolution. Simulations including ram pressure 
show oscillations in their SFR, which can be traced to density 
fluctuations in the galactic disk. A possible cause is the numeri¬ 
cal problem in SPH (see Sec. 3.3), seeing the compressions also 
occur when the star formation and galactic-wind subgrid model 
is disabled. The oscillations, however, are not comparable to the 
physical effects that can be seen in our plots. 

In Fig. 4, the second and highest peak at a SFR of 19 M Q yr _1 
coincides with the time of the galaxies’ first encounter (t ~ 
0.4 Gyr). It is present in both curves that show the merger SFR 
in the ICM and the vacuum merger’s SFR. As mentioned in the 
introduction, the origin of these peaks is mostly the triggered 
gas inflow to the center and higher gas densities due to the dis¬ 
turbance induced by the merging partner. 

To put the mergers into perspective, the added vacuum star 
formation of the individual galaxies involved in the merger (i.e., 
added SFR of the two separate simulations p 1 and p6) is also 
shown in Fig. 4. In a vacuum, the non-merging galaxies form 


stars at a slowly declining, combined rate of around 3 Moyr^ 1 . 
The vacuum merger has comparable SFRs up until the first en¬ 
counter of the merging galaxies at t — 0.4 Gyr and then peaks 
sharply at 13 M G yr~’. This increase surpasses the combined rate 
of the individual galaxies in the ICM without a merging part¬ 
ner (sum of SFRs for the simulations labeled “in 1” and “in 6”), 
which is the top curve in Fig. 4. The first encounter peak in 
the ICM is 25% higher than in the vacuum merger and higher 
than the non-merging galaxies in ICM. This collaborative ef¬ 
fect is a result of the core gas being more compact under ram 
pressure, making it more susceptible to gas inflows triggered by 
the merger. The effect, however, is short-lived: The second en¬ 
counter peak at t « 0.86 Gyr is small in comparison to the vac¬ 
uum case and barely visible in the ICM merger. By the time of 
the second encounter, >25% of the gas is stripped, and the gas 
and stellar disks of both galaxies are severely disturbed by the 
interaction. These results are consistent with the results of K8. 
We use more moderate parameters for the galaxies’ gas fraction 
(15% instead of 25% in K8), resulting in an overall lower SFR 
compared to K8 and K9. For otherwise comparable galaxies, 
they have peak values of 26 M G yr 1 at an enhancement in ICM 
of a factor of 3.5 to 8. The standard case represents the highest 
relative peak enhancement of SFR in our simulations. 

The impact of changes in the ICM density is shown in Fig. 
5. We omit the plots for wp21 and wp29, i.e. highest and low¬ 
est simulated ICM densities. They do not differ significantly 
from their comparison runs for “ram pressure” - only (ip\p21 + 
ip6p21 in case of wp21) or vacuum (mp6) in the case of wp29 - 
and thus show a limit for possible collaborative effects. 

In both the ICM cases in Fig. 5a (merging and non-merging), 
overall SFR is steeply declining for picm = 5 x 10 2S g cm~ 3 . By 
the time of the first encounter, when the first SFR peak appears in 
the vacuum merger, more than 50% of the gas is stripped. Thus, 
we barely see an enhancement in the ICM merger’s peak SFR 
over the non-merging galaxies in ICM in Fig. 5a. The vacuum 
merger peak in SFR at the first encounter in Fig. 5a is below 
the value for the non-merging galaxies in ICM. In contrast, the 
lower ram pressure runs, wp529 and wstd, show a distinct first 
encounter peak in the ICM merger. Indeed, it is a handy heuristic 
for the simulations in this work: if any encounter peak in the 
vacuum merger is at least comparable in height to the overall 
SFR in ICM single galaxies, the peak can also be seen in the 
ICM merger. 

Varying the relative speed between ICM and the galaxies 
(iwr500) showed only a decrease in overall star formation by 
about 15 %, when compared to the standard case. Otherwise, 
vvit500 showed the same peak heights and general behavior as 
the standard case and is therefore not shown in detail. This be¬ 
havior is expected, e.g. from K9, because the ICM density is the 
more crucial parameter for SFR in traditional SPH simulations 
than is the speed relative to the ICM, although the lower ram 
pressure via v re i does affect the stripping rate (see Sec. 3.2). 

The SFRs for different mass ratios are shown in Fig. 6. 
Since the partner is more massive and can trigger a stronger star- 
burst, the simulation with mass ratio 1:4 (Fig. 6a) has more pro¬ 
nounced merger peaks in vacuum than the reference case. In 
the ICM merger, this translates into slightly higher merger peaks 
as well, but not a proportionate increase. With p = \ : 8 in 
Fig. 6b, we seem to be at the limit mass ratio for the merger to 
show any changes in SFR in ICM environments. In the counter¬ 
rotating runs, such SFR plots proved to be very similar to their 
co-rotating counterparts in Fig. 6, with the small increase in 
SF expected from theory (stemming from additional instability). 
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t [Gyr] 

(a) ICM density p IC M = 5 x 1(T 28 g cnr 3 



(b) ICM density p IC m = 5 x 10 29 g cm 3 


Fig. 5: SFRs for different ICM densities. The higher density 
run wp528 in (a) shows a minute merger peak and is otherwise 
indistinguishable from the single galaxies in ICM. Run wp529 in 
Fig. (b), on the other hand, has a visible first encounter peak in 
ICM that is slightly smaller than in wstd but that still surpasses 
the peak value in a vacuum. 


They were omitted for brevity, but included in the summarizing 
Fig. 8. 

Other configurations of angular parameters, i.e., ft, y (see 
Fig. 3), and the ICM direction are not all shown in distinct plots, 
but are part of Fig. 8. Variations in these parameters flatten the 
merger-peaks in vacuum, because the interaction time between 
the two disk increases or the area for ram pressure to act on is 
reduced. Different displacements (impact parameters) and initial 
merger velocities, b and v mer , show little or no merger peaks in 
the ICM, since the interaction is too short. An example of this 
general trend is shown in Fig. 7 for y = 90°. Here the smaller 
disk feels an edge-on wind, with the larger encountering the ICM 
face-on. The vacuum merger peaks at around t = 0.4 Gyr and 
1 = 0.85 Gyr are broader compared to wstd but less pronounced. 
In ICM they are broader as well but barely stand out from the 
general enhancement due to ram pressure. Of these different 
geometric configurations, only wbl 0 has a sharp first encounter 
peak at 17 M Q yr _1 (15 M Q yr 1 in the vacuum merger). Just 
like wstd, the smaller galaxy hits the larger relatively close to 



(a) Mass ratio p = 1:4 



(b) Mass ratio p = 1:8 

Fig. 6: SFRs for different mass ratios between the merging 
galaxies. Fig. (a), showing wp4, has a more pronounced sec¬ 
ond encounter peak in both ICM and vacuum compared to wstd 
and more even so to wp8 in (b). 



Fig. 7: Simulation ivy 90 and related SFRs with inclination angle 
of the smaller galaxy y = 90°. The galactic disks are oriented 
perpendicular to each other, resulting in a longer crossing time 
and thus broader but flatter vacuum merger peaks. 
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Fig. 8: Total mass of newly formed stars relative to the added total new star mass of the merger constituents in vacuum. Values are 
at t — 1 Gyr (end of simulation time). Mergers in the ICM show slightly but not significantly more star mass than the individual 
galaxies. 


its center in a face-on collision and thus yields a comparable 
outcome. 

Finally, Fig. 8 shows the relative mass of new stars at the end 
of the simulations relative to the corresponding mass in the in¬ 
dividual galaxies in vacuum. Only one case (wvm!5 with lower 
relative merger velocity) has formed noticeably fewer stars in the 
ICM merger than its constituent non-merging galaxies in ICM. 
We attribute this slightly higher SF to stronger perturbations and 
triggered inflows to the star-forming center due to the merger, 
as well as the potential wells being deeper (added mass com¬ 
pared to single galaxies). The variance in overall star formation 
with the individual merger and ICM parameters follows no clear 
trend and is comparably small. It ranges between almost no dif¬ 
ference in the case of y - 45° to a difference in star mass of 
18 % between ICM single galaxies and ICM merger in the sim¬ 
ulation with edge-on merger and edge-on ICM (wee). The total 
mass fraction of new stars in the ICM mergers does not follow 
vacuum merger tendencies; e.g., relative mass of new stars in¬ 
creases slightly from mvmSO to mvml5 (vacuum mergers) but 
decreases between the ICM case (wvmSQ, wvmlS). 

With the exception of a narrow time frame around the first 
impact of the two merging galaxies, our simulations show neg¬ 
ligible contributions of minor mergers to the SFR even in low- 
density environments (e.g., on the outskirts of typical clusters, 
outside the virial radius). 

3.2. Stripping and star formation in the wake 

The time evolution of gas and star-mass fractions in the wake is 
plotted in Fig. 9 for the reference case wstd. We show these 
same fractions for the combined masses of the non-merging 
galaxies as well, in order to provide a comparison. With the 
star formation recipe used, wake stars form mainly in nodes of 
stripped gas (see SI, K8). Our criterion for wake stars and gas is 
simple: We consider particles to be part of the wake, if they lag 
40 kpc behind the common center of mass of the stellar disks. 
This way we can use a uniform criterion across our simulations, 



Fig. 9: Fraction of stripped gas to total gas mass of the galaxies 
(upper two lines) and mass of stars newly formed in the wake 
relative to the total mass of new stars (lower two lines). Shown 
are the values for both the merger (solid lines) and the added 
single galaxies (dotted lines) in ICM for wstd. 


thereby covering all the unbound gas. The result of a binding 
energy-based criterion for identifying wake gas agrees reason¬ 
ably with this simpler criterion, with a delay of ^50 Myr. 

Initially all curves in Fig. 9 show a fast increase in wake 
mass fractions. When the gas stripping approaches the equi¬ 
librium predicted by Gunn & Gott (1972), the slope becomes 
shallower. In Fig. 9 this occurs after 1 ~ 0.35 Gyr for gas and 
1 ~ 0.5 Gyr for the stellar component. In principle, the stellar 
component follows the trend set by the galaxies’ gas with a tem¬ 
poral delay of just over 0.1 Gyr. The increasing offset is a result 
of the reaccretion of wake stars by the galaxies, because stars 
no longer feel the ram pressure caused by the relative motion 
against the ICM. Notable as well is a difference in wake mass of 
both the stellar and gaseous components, when comparing single 
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Fig. 10: Stripping fraction at a merger mass ratio of p — 1 : 8. 
The curves shown here for the simulations w//8 in Table 3 are 
quantitatively the same as in Fig. 9 with slightly more gas being 
stripped. 


and merging galaxies. After the first 0.3 Gyr, the single galax¬ 
ies experience a higher stripping fraction than in the respective 
merger scenario. The added mass of the second galaxy in the 
merger, as well as the fact that the merging partner collects some 
of the stripped gas after first passage, explains this discrepancy. 
The additional galaxy increases the potential barrier for gas to 
escape, once the galaxies are sufficiently close. The companion 
galaxy overtakes the stripped gas in many of our simulations and 
binds some of it again. 

In our mass regime and predominant merger geometry, this 
effect compensates for more turbulent gas and tidal features of 
the merger (such as the bridge), which are more readily stripped. 
It does, however, depend on the exact merger geometry, as can 
be seen from wee and wfiy, where re-accreation of gas hardly 
occurs. Figure 10 shows the same plot for a smaller secondary 
galaxy. It is comparable to that of Fig. 9 with a slightly higher 
stripping fraction and higher ratio of star formation in the wake, 
although it shows a jump in stripping fraction for single galaxies 
at ts;0.65 Gyr. The cause of this jump is the core of the smaller 
galaxy finally being stripped. It is only weakly bound until that 
time. In the ICM merger, the larger galaxy helps retain this core. 

Figure 11 shows the situation at higher density (Fig. 11a) and 
lower density (Fig lib). In case of higher density (Fig. 11a), 
the stripping rate is enhanced by about a factor of three com¬ 
pared to wstd. About 40% of stars are located in the wake by the 
time saturation is reached in the higher density ICM. The rela¬ 
tive difference between the curves for single and merging galax¬ 
ies becomes smaller at higher densities. At higher densities, the 
stripping is faster and more severe, giving the merger less time 
to influence the process. Clearly, ram pressure is the dominant 
effect here. 

At lower densities (Fig. 1 lb), gas stripping is a more con¬ 
tinuous process, with no clear saturation phase. The instability 
caused by the merger enables gas to be stripped even at later 
times, when the single galaxies would have reached their satu¬ 
ration point. In the merger, however, gas stripped at later times 
forms fewer stars than do single galaxies. The disturbance in 
the gas disk is much greater then, allowing for less star-forming 
nodes in the wake. At even lower densities (simulation wp29; 
plot omitted), only 5% of gas is stripped by the end of the sim- 



t [Gyr] 

(a) ICM density p IC M = 5 x 10 -28 gem -3 



Fig. 11: Stripping fractions for the different ICM densities of 
the simulations wp528 and wp529. Note the different scales on 
the ordinate. The case for lower density (b) does not reach its 
saturation point within the simulation time of 1 Gyr and has a 
later onset in the wake star fraction. 


ulation time with less than 1 % of all newly formed stars in the 
wake. 

In the graph of Fig. 12, where the stripping fractions of the 
runs with relative velocity of v r = 500 l< m s 1 are shown, we 
see drastically decreased stripping and even less star mass in the 
wake. Contrary to the situation for overall SFR in the previous 
section, a lower velocity relative to the ICM has a strong influ¬ 
ence on the location of star formation (wake vs. disk). Toward 
the end of the simulation in Fig. 12, after t = 0.75 Gyr, even 
re-accretion of stars is visible, as the lower solid curve declines 
slightly. This is the only simulation to show a decrease in the 
wake’s stellar mass fraction. 

Figure 13 concludes this section by summarizing the results 
of all simulations. In Fig. 13a one can see the amount of gas 
in the wake, again relative to the respective galaxies’ total gas 
mass. Shown are only the values at the end of the simulation. 
Geometric merger parameters show little effect on the amount of 
gaseous material in the wake between merging and non-merging 
galaxies. Compared to the single galaxies in ICM, slightly less 
material gets stripped in the scenarios involving a minor merger, 
though this is dependent on merger geometry. The difference in 
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(a) Gas fraction in the wake 


(b) Fraction of new stars in the wake 


Fig. 13: Gas and new star mass in the wake relative to the respective components’ total mass (t — 1 Gyr). Shown are the values for 
both the merger simulations in ICM and the added single galaxies in ICM that make up the respective merger. 



Fig. 12: Stripping fraction for a lower relative velocity of v r = 
500 kms^ 1 than in wstd of Fig. 9. The plot shows that sig¬ 
nificantly less mass is being stripped, and stars are noticeably 
re-accreted toward the end of simulation. 


stripped mass fraction between merging and non-merging galax¬ 
ies in our simulations ranges from almost none (w[3y) up to 25%; 
for example, w//8 has 26% of its gas residing in the wake, while 
its constituent galaxies i\i 8 and ifj.1 combined have 33%. 

Figure 13b shows the same situation for the fraction of new 
stellar mass in the wake to total new stellar mass by the end of 
the simulations. Overall, the fraction of new stars formed in the 
wake is lower in the merger simulations compared to the single 
galaxies. As little as half the fractional new stellar content in the 
wake of single galaxies can be seen in the wake of the respective 
merger cases. This is to some extent a result of the added mass 
in the mergers. More so, in most of our simulations the merging 
partner comes close to the wake of the other galaxy. That en¬ 
ables the partner to accrete some of the stars from the wake of 
the other galaxy. In the runs where the merging partner does not 
approach the wake (wee and wnifiy) closely, the differences are 
smaller. At higher densities, stripping is mostly complete before 
merger effects start to contribute. At lower densities, ram pres¬ 
sure stripping is no longer a factor, and the difference between 


merging and non-merging galaxies is small. Only intermediate 
densities and favorable trajectories of the merging partner lead to 
collaborative effects. In the light of our simulations, taking mi¬ 
nor mergers into account in the process of ram pressure stripping 
does not yield significantly different results from ram-pressure- 
only scenarios (except for the wake’s stellar content in the men¬ 
tioned circumstances). 

3.3. Discussion and related work 

Unlike our simulations of constant ram pressure, real cluster 
galaxies experience varying picm, depending on their path within 
the cluster. The recent study of Bekki (2014) uses a /j-pro file for 
Picm and model ram pressure along different galaxy paths with 
an SPH code to quantify the effect of varying density on SFRs. 
Their “MW-type” galaxy model is similar to the larger merg¬ 
ing partner in our simulations with less total mass (10 12 M 0 ). 
Their “M33” is most similar to our ip 8. The Virgo and group¬ 
like ICM models in Bekki (2014) yield a SFR enhancement by 
a factor of 2-3 (compared to our 3-5x enhancement), but only 
in the Gyr after the pericenter passage of the galaxies and very 
shortly before. At this time in the galaxies’ orbits within the 
Virgo model, ICM density would be > 5 x 10~ 28 g cm -3 . This is 
different behavior from our individual galaxies in ICM that reach 
their enhancement sooner, after only 0.2 Gyr and generally show 
stronger enhancement. 

A physical explanation may be that gradual infall removes 
gas gently without causing the conditions that increase SFR in 
our simulations. With star-forming gas then missing at later 
times, enhancement in SFR at higher density might not be as 
drastic until the short-term peak density at pericenter passage 
comes into play. It is difficult, though, to separate the physi¬ 
cal effects from the different numerical models for star forma¬ 
tion used in both our studies and other details of the simulation. 
Jachym et al. (2007, 2009) studied stripping rates using the same 
approach for varying picm as Bekki (2014). Their results indi¬ 
cate that short time variations in ram pressure do not affect the 
stripping rate. Instead, the total encountered column density of 
ICM is the crucial parameter for the amount of material stripped. 

We identified the deeper potential well and re-accretion by 
the merging partner’s increased reach to be the cause of less 
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pronounced gas stripping and fewer stars in the wake of galax¬ 
ies merging in the ICM. We expect qualitatively similar strip¬ 
ping fractions under variable ram pressure. When compared to 
simulations with adaptive mesh refinement (AMR), our simu¬ 
lations, and SPH simulations in general, differ greatly in SFR. 
AMR simulations consistently show quenching rather than an 
increase in SFR under ram pressure. Tonnesen & Bryan (2012) 
compare the approach of K9 to their AMR simulations and to 
observations. They attribute SPH’ problems with resolving fluid 
instabilities and the resulting under-mixing, as well as the stiff 
equation of state as a cause for a high SFR in the tail. These 
and related problems with the standard pressure-entropy formu¬ 
lation of SPH have been discussed in Agertz et al. (2007) and 
Hopkins (2013). The problems are most notable in the stripped 
tail, where much of the star formation takes place in dense nodes 
of gas. This leads to increased stability and higher density of 
the star-forming nodes in the wake. Generally in GADGET-2’s 
formulation of SPH, strong density discontinuities lead to a false 
pressure estimate (which relies on a smooth density gradient), re¬ 
sulting in an unphysical separation force. HeB & Springel (2012) 
showed that a gap between ICM and ISM can form because of 
that, which could cause the density oscillations in the disk. This 
also keeps disk particles in the star-forming regime, especially at 
later times, when fluid instabilities and mixing would decrease 
star formation. Observations support lower SFRs than in our 
simulations (Boissier et al. 2012). 

Vollmer et al. (2012) finds possible enhancement factors of 
1.5 - 2 on the windward side of his sample galaxies in Virgo, 
while the overall effect is suppression of star formation, espe¬ 
cially in the tail. Star-forming tails were found, among others, 
by Sun et al. (2007) and more recently Ebeling et al. (2014). 
The former estimates the total mass of current star bursts in the 
tail’s nodes to be several 10 7 M 0 , asserting the total stellar mass 
formed in the tail could be a few times higher from already faded 
nodes. Ebeling et al. (2014) shows strong galactic starbursts that 
are likely to be very short in duration. Hester et al. (2010) inves¬ 
tigated the Virgo dwarf galaxy IC3418, finding nodes of dense 
gas, generally rare in other Virgo galaxies, which Jachym et al. 
(2013) studied in greater detail. The total stellar mass in the 
tail of IC3418, however, is less than 1% of the galaxies stel¬ 
lar mass with tail SFRs of ~ 0.002 M Q yr _1 . HeB & Springel 
(2012) compared the behavior of galaxies in a windtunnel be¬ 
tween GADGET-2 and the moving mesh code AREPO (Springel 
2010). They demonstrated an overestimation by a factor of two 
in SFR and a smaller stripping fraction in SPH. 

Our objective with this paper has been to determine whether mi¬ 
nor mergers have an effect on stripping and SF that should be 
considered. Although the issues discussed with SPH pose con- 
ceptional problems, we use relative quantities and compare, for 
instance, the relative increase in stripping fraction or SF between 
merging and non-merging galaxies in the ICM. Moreover, our 
conclusion is that the effects of merging in the ICM are not sig¬ 
nificantly different from regular ram-pressure effects on single 
galaxies. 

Hayward et al. (2013) found that pure merger simulations 
using the standard density-entropy formulation of SPH com¬ 
pare well with results from state-of-the-art, grid-based hydro- 
dynamic codes (i.e., AREPO). The physical processes that could 
conceivably have led to synergies of ram pressure and merging 
(gas inflow to the galactic center, shocks, and disturbance of the 
gaseous disk from the merger, augmented by increased pressure 
from the ICM) are accurately modeled in the case of pure merg¬ 
ers or overestimated in SPH or while the counteracting stripping 
and mixing of the gas phases is underestimated. Therefore the, in 


any case, small collaborative effects should be even less distinct 
in these codes, and our conclusion should hold. 

4. Conclusions 

We have presented GADGET-2 simulations of merging galaxies 
in the ICM and compared them to galaxies evolved without a 
merger partner. Our simulations involved gas-rich spiral type 
galaxies in minor mergers (mass ratios between 1:4 and 1:8) and 
included their first two encounters. 

In general, the merging galaxies retain 5% to 10 % more gas 
in their disk than do galaxies without a merger partner; i.e., less 
of their gas is located in the ram-pressure stripped wake. The 
added mass of the merging parter and re-accretion of wake gas 
by the partner (depending on trajectory) compensate for the dis¬ 
turbance of the gas disk that should make stripping easier. Ow¬ 
ing to re-accretion of stars and the added gravitational potential 
of the merging partners, we typically found 50% fewer stars in 
the wake of merging galaxies. This effect depends strongly on 
the trajectory of the galaxies with respect to the wind direction. 
The relative differences are smaller in simulations with higher 
ICM density as ram pressure begins to dominate. 

Even on the outskirts of clusters, with typical densities of 
Picm = 10 28 g cm -3 , our simulations indicate little overall ef¬ 
fects on SFR owing to ongoing merging activity. We observe 
a short-term enhancement in SFR, compared to the merger in 
vacuum of up to 50 %, but typically only 20 %. These enhance¬ 
ments occur at the time of the first fly-through, when ram pres¬ 
sure condensed the core gas, making it more susceptible to SF 
induced by the merging partner. Second encounters rarely show 
any SFR peaks in most cases. Where they do, the enhancements 
are not higher than 20 % over the added SFR of the individual 
galaxies in ICM. By the time of the second encounter, the galax¬ 
ies’ disks are too disturbed, and too much gas is stripped. In any 
case the short-term enhancements do not last longer than 50 Myr 
and do not contribute much to the overall star mass. When com¬ 
paring the new star mass formed by the end of our simulations, 
generally the merging galaxies have formed slightly more stars, 
typically around 5 % but up to 18 % more than the individual 
constituent galaxies in ICM after a time of 1 Gyr. Depending 
mostly on ICM density and geometric parameters, the reason for 
this higher star formation in our simulations is a denser core and 
gas inflow triggered by the merging partner. We discussed possi¬ 
ble issues in Sect. 3.3. They boil down to unphysical pressure es¬ 
timates at the interface between ICM and ISM, possibly leading 
to an overestimation of the collaborative effects. For overall star 
formation, the effect of varying geometric and ICM parameters 
related to ram pressure stripping dominates those of the minor 
merger parameters (mass ratio, merger geometry), and collabo¬ 
rative effects are small in any case. 
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